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We argue that thermal machines can be understood from the perspective of Virtual qubits' at Vir- 
tual temperatures': The relevant way to view the two heat baths which drive a thermal machine is as 
a composite system. Virtual qubits are two-level subsystems of this composite, and their virtual tem- 
peratures can take on any value, positive or negative. Thermal machines act upon an external system 
by placing it in thermal contact with a well-selected range of virtual qubits and temperatures. We 
demonstrate these claims by studying the smallest thermal machines. We show further that this per- 
spective provides a powerful way to view thermodynamics, by analysing a number of phenomena. 
This includes approaching Carnot efficiency (where we find that all machines do so essentially by 
becoming equivalent to the smallest thermal machines), entropy production in irreversible machines, 
and a way to view work in terms of negative temperature and population inversion. Moreover we in- 
troduce the idea of ''genuine'' thermal machines and are led to considering the concept of "strength" 
of work. 



I. INTRODUCTION 

The field of quantum thermodynamics fP-Sl has made 
great strides in gaining a basic understanding of phe- 
nomena at the intersection of quantum mechanics and 
thermodynamics. This line of research goes all the way 
back to the 1950's w^hen thermodynamic analysis of 
lasers was investigated [4-6J. Since then there has been 
significant interest in many related areas, w^ith many re- 
sults on quantum thermal machines [7-16J, finite-time 
thermodynamics |[T7ti22l and the second law^ li23rr26J 
(and references therein). 

Of particular interest to this w^ork is one recent devel- 
opment, namely the study of the smallest possible self- 
contained thermal machines [27 -31 J. By self-contained 
we mean that no sources of external work or other form 
of control are allowed; only incoherent interactions with 
thermal baths at various temperatures. These smallest 
self-contained thermal machines are arguably the most 
elementary thermal machines. As such, they are inher- 
ently simple and transparent, they offer us a view into 
the core of thermodynamics unobstructed by unneces- 
sary details and complications. What we want to show 
here is that because of this, their study leads to a new 
view of thermodynamics, allowing for general conclu- 
sions to be drawn about the way thermal machines ulti- 
mately function. 

At the core of this view of thermodynamics lies the 
notion of a virtual qubit and its virtual temperature. A vir- 
tual qubit is a two-level subsystem of the two baths that 
drive a thermal machine, when considered as a compos- 
ite system. Different virtual qubits have different vir- 
tual temperatures. Using these concepts, we will show 
first that the smallest thermal machines - refrigerators, 
heat pumps and heat engines - all function via a sim- 
ple mechanism; they place an external system in thermal 
contact with a well chosen virtual qubit, at a well cho- 



sen virtual temperature. More complicated thermal ma- 
chines, including classical ones, use essentially the same 
mechanism. The only difference is that they couple the 
external system to many virtual qubits, at many virtual 
temperatures, all at the same time. 

We will demonstrate further the breadth of this new 
notion. First of all the virtual qubit provides a natu- 
ral way to understand work, in terms of population in- 
version and negative temperatures. This in turn shows 
that work has an additional property, which we term 
"strength''. Additionally this notion is seen to be partic- 
ularly powerful in understanding efficiency, especially 
the Carnot limit. The smallest thermal machines will be 
shown to have an efficiency which is always universal 
(independent of model details), whilst any thermal ma- 
chine which approaches the Carnot efficiency functions, 
essentially, as the smallest machines in the limit, util- 
ising only a single virtual temperature. Moreover the 
strength of work will be seen to vanish in the Carnot 
limit, reinforcing the weak nature of the Carnot limit. 
Finally, we will see that the resources a pair of thermal 
baths provide is captured by the notion of virtual tem- 
perature, which in turn allows for the idea of a 'genuine 
thermal machine' to be introduced. 

The paper is organised as follows: In section |ll] we 
introduce the virtual qubit and its virtual temperature 
and outline in section |lll] how this notion is used to un- 



derstand thermal machines. In section IV we begin to 
demonstrate explicitly this idea for the smallest thermal 
machines, discussing first refrigerators and heat pumps 
in section lyj and then work and heat engines in section 
^ In section |VII| we introduce the 'strength' of work 



VI 



and its manifestation for both finite-dimensional and 
infinite-dimensional systems. We begin our discussion 
of efficiency in section VIII[ where we discuss the uni- 
versality of the efficiency of the smallest machines. Con- 
tinuing in section [lX| we discuss how thermal machines 
approach the Carnot limit, whilst in section |X| we study 
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FIG. 1: Schematic diagram of bath virtual qubits. 

(a) Two baths, one at temperature Ti and one at tempera- 
ture T2. In each bath we identify two arbitrary energy eigen- 
states, and with spacing Sj. (b) These two baths 
can be viewed as a single composite bath. Two energy eigen- 
states of this composite system are |G)g^ = \^)b^\s)b2 

= \g\^ . The Virtual temperature' % of this Virtual 
qubit' can take on any value - positive or negative - depend- 
ing on Si and 82- 



entropy production in irreversible machines. Finally in 
section XI we introduce 'genuine' machines. 



11. THE VIRTUAL QUBIT 

Suppose we have tw^o non-interacting thermal baths 
at different temperatures, Ti and T2. The state of each 
bath by itself is rather trivial, being simply a Boltzman- 
nian distribution. However, thermal machines work by 
accessing both baths. It is therefore far more relevant to 
look at the two baths together, as a composite system. 
The structure of this system, which may seem no more 
complicated than its constituent parts, has in fact a rich 
structure. A 'Virtual qubit" is the most elementary sub- 
space in the Hilbert space of the composite bath, i.e. a 
two-dimensional subspace. Most important are virtual 
qubits in which the two (basis) states are energy eigen- 
states of the composite system. What makes the virtual 
qubit interesting is that in general it behaves as if it has a 
virtual temperature which can be very different from 
either bath temperature, Ti or T2; in particular % can be 
much smaller, or larger than either, or even negative. 

To be more precise, let us consider two arbitrary en- 
ergy eigenstates in bath 1 with energy spacing <f 1, which 
we denote \g)^^ and \e)^^. Similarly, consider two states 
in bath 2 with energy spacing 82^ denoted \g)^^ and 
\e)^^. These are depicted in FIG. [l] (a). From these 
states we can identify four energy eigenstates of the 
composite system, namely Ig)^^, Ig)^^, Ig)^^ 
and Of particular interest are virtual qubits 

formed of the states \e)g^ |g)g^ and \g)^^ Without loss 
of generality we may take S2 > £1. With this choice, 
these two states form a virtual qubit with energy spac- 
ing = £2 — £i- It is convenient to define the virtual 



ground and excited states as, respectively 

depicted in FIG. IT] (b). The "virtual temperature" of this 
qubit can be found by looking at the ratio of populations 
and p^, of the excited and ground state. 
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Since each bath is individually in a Boltzmannian distri- 
bution, we know that the populations of the eigenstates 
from a single bath satisfy the relations 



Therefore we see that 
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as the virtual temperature of the virtual qubit. There are 
a number of important observations which we can now 
make. 

The first is that this virtual temperature can take a 
range of values, crucially all temperatures outside the 
range of Ti and T2 - it can be smaller or larger than 
both, and can even take negative values. (Note that if 
we had instead considered virtual qubits formed of the 
states |g)g^ \s)b2 I^)bi I^)b2 then the virtual tempera- 
tures would have been limited to lie between Ti and T2.) 

Second is that in a pair of macroscopic baths there will 
be essentially infinitely many virtual qubits at any given 
spacing S^, and any given temperature It is precisely 
this fact that gives the composite system its thermody- 
namic power. 



III. THERMAL MACHINES 

The central idea of this paper is that: 

All a thermal machine does is to place an external 
system in direct thermal contact with a restricted 
set of virtual qubits in the composite hath, having 
only a restricted range of virtual temperatures. 
The external system simply reacts to the virtual 
qubits as it would react if put in thermal contact 
with real qubits having the same virtual temper- 
atures. 
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If the temperatures are predominantly smaller than 
those of the individual baths, the machine is a refriger- 
ator. If the temperatures are predominantly larger, then 
it is a heat pump. Finally, if the temperature is nega- 
tive, then the machine is a heat engine. In this sense 
thermal machines act simultaneously as a coupler and 
a filter - they provide a coupling between the system 
and the thermal reservoirs, but also filter out only a re- 
stricted range of virtual temperatures. 

The important question is how exactly does a ther- 
mal machine achieve this task. In particular, how can 
a thermal machine filter out selected virtual qubits from 
the composite bath? This ability seems surprising since 
we can construct thermal machines without detailed 
knowledge of the spectrum of the bath. In the next sec- 
tion we will give an answer to this question for the case 
of the smallest thermal machines. In fact, as we will see, 
the smallest thermal machines couple not to a range of 
virtual temperatures, but to only a single one. In this 
sense, they act as perfect filters and are therefore the 
"cleanest'' thermal machines. We will see later in Section 
IX that the ability to perfectly filter is intimately related 
to the ability to reach Carnot efficiency. 



IV. HOW THERMAL MACHINES WORK: THE 
SMALLEST THERMAL MACHINES 

A. The machine virtual qubit 

All of the smallest thermal machines have at their core 
the same basic mechanism. Each machine consists of 
two qubits, qubit 1 that is in thermal contact with bath 1, 
and qubit 2 that is in thermal contact with bath 2. These 
two qubits will then interact with each other and an ex- 
ternal system via an interaction Hamiltonian. It is illu- 
minating to consider the state of the two qubits before 
this interaction is turned on. Suppose that qubit 1 has 
energy spacing Ei and qubit 2 energy spacing £2- The 
explicit choice of these energies is part of the design of 
the machine and is chosen by us. In the absence of in- 
teraction with the external system, each qubit interacts 
only with its own thermal bath. As such, each will reach 
thermal equilibrium at the corresponding bath temper- 
ature. 

For simplicity, throughout this paper we will only 
consider the case of weak coupling between the machine 
and the bath, as this will allow us to focus on resonant 
(i.e. energy conserving) interactions only. In this regime 
the thermal state of each qubit will be a Boltzmannian 
distribution with the Hamiltonian being the free Hamil- 
tonian of the qubit. Thus the thermal state of each qubit 
is given by 
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where Hj = E/|l).(l| is the free Hamiltonian of each 
qubit. 

Let us look at this trivial thermalisation process, in 
which each qubit simply reaches equilibrium with its 
own bath, from a different, and ultimately more rele- 
vant, point of view. The two qubits have 4 energy eigen- 
states, |0)J0)2, |0)Jl)2, |1)J0)2 and |l)Jl)2. Of particular 
interest are the two states |0)Jl)2 and 1 1)^0)2. In sim- 
ilar fashion to the previous section, let us interpret the 
Hilbert space spanned by these two states as a virtual 
qubit, with ground and excited state 



|0)„ = |1)J0), 
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This will be our machine virtual qubit. The energy spacing 
of these two states is E-^ = E2 = Ei. By repeating a 
similar analysis to that carried out for the virtual qubit 
in the composite bath, we find that this machine virtual 
qubit reaches the temperature Ty 
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We can interpret this result as saying that the machine 
virtual qubit entered into thermal contact with virtual qubits 
in the composite bath and thermalised to their temperature. 

Seen from this point of view, the thermalisation pro- 
cess may seem surprising. As we saw previously, there 
are many virtual qubits in the composite bath having the 
same energy-level spacing as the machine virtual qubit, 
i.e. satisfying E^ — — — It would appear there- 
fore that all such virtual qubits should couple resonantly 
to the machine virtual qubit. However, the temperature 
of the virtual qubits depends not only on the level spac- 
ing E^ — E\, but also on the separate values of and E\) 
coupling all of these to the machine virtual qubit will re- 
sult in it attaining some average temperature. The key 
to understanding the behaviour of the machine virtual 
qubit is to realise that the interaction with the composite 
bath is not arbitrary; it occurs only via the local couplings 
of qubit 1 to bath 1 and qubit 2 to bath 2. Each machine 
qubit will only interact resonantly with a corresponding 
qubit in the bath, so E\ = E\ and E^ = £2- This allows 
the machine to couple its virtual qubit in a very selec- 
tive way to the virtual qubits in the composite bath and 
to single out the precise temperature T^. In the termi- 
nology of the previous section, we see that it is the local 
resonant couplings which provide the filtering for these 
machines. Furthermore, we see that only a single is 
selected, and as such this machine acts as a perfect fil- 
ter. Finally, we note that the specific value of filtered 
out is chosen at our disposal by the the way in which we 
engineer the machine, i.e. by the choice of E\ and £2- 



B. Coupling to the external system 

We now have all the necessary tools to understand the 
behaviour of these thermal machines: All we have to do 
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FIG. 2: Schematic diagram of the smallest thermal machines. 

(a) The smallest machine comprises two qubits, one with en- 
ergy spacing E\ in thermal contact (wavy line) with a bath at 
temperature T\ and the other with energy spacing £2 iri ther- 
mal contact with a bath at temperature T2. (b) The ultimately 
more relevant way to view this system is that the 'machine 
virtual qubit' with spacing = £2 — £^1/ is in thermal contact 
with virtual qubits in the composite bath at the virtual temper- 
ature Ty. The local energy-level spacings £1 and £2 allow the 
machine virtual qubit to 'filter' out a single virtual tempera- 
ture from the composite bath. 



now is to place the external system in thermal contact with 
the virtual qubit of the machine. 

As an example, let us consider the simplest case, 
w^here the system w^hich the machine will act upon is 
itself another qubit. This qubit shall be called qubit 
3 and shall have energy spacing £3. Given this spac- 
ing, w^e ensure that the machine is engineered such that 
£2 — £1 = £3- With this choice, we see that the machine 
virtual qubit and the external system have equal energy- 
level spacings. It is therefore possible to introduce an ar- 
bitrarily weak interaction which allows the virtual qubit 
and external system to resonantly exchange energy - i.e. 
we place them in thermal contact. 

This can be done via the interaction Hamiltonian 

H/„f =g(|0)Jl)JO)3/l|,(0|3(l| + |1)J0)J1)3/0|,(1|3(0|) 

(10) 

which was previously introduced in [27J and in- 
duces transitions between the two degenerate states 
1 0) J 1)2 10)3 ^ 1 1) J 0)2 11)3. Again the most relevant way 
to write this Hamiltonian is in terms of the virtual qubit. 



Hint = g{\Ol\lUMP\ + |l)JO)3iOL,(l|) (11) 



which induces transitions between the two degenerate 
states 10)^11)3 ^ 11)^10)3. In this notation is becomes 
manifest that the interaction Hamiltonian is the one 
which generates the unitary SWAP between the external 
system and the machine virtual qubit. Hence this inter- 
action places these two qubits into thermal contact. 



C. Summary 

We now reach the crux of our analysis. We see that the 
way in which the whole process works is that selected 
virtual qubits of the composite bath, at a selected virtual 
temperature, are placed in thermal contact with a virtual 
qubit of the machine. In turn this virtual qubit of the 
machine is placed in thermal contact with the external 
system. Hence all that happens is an ordinary thermal- 
isation process of the external system with the virtual 
qubits in the composite bath, mediated by the machine 
virtual qubit. If the external system is otherwise isolated 
from an external environment, then it reaches the tem- 
perature (for an explicit calculation see Appendix [A} 
and FIG.|3|. On the other hand, if the system is also in 
contact with an external environment, then it will reach 
some (stationary) state similarly to any system in con- 
tact with two thermal baths, see FIG.|4] 

We conjecture that this process is in fact a special case 
of the general procedure which all thermal machines 
use. However, a thermal machine generally has many 
different energy levels, hence many different machine 
virtual qubits. In turn they couple to many different vir- 
tual qubits in the composite bath, covering a range of 
virtual temperatures, instead of a single one, thus com- 
plicating the situation. It is also further complicated by 
the fact that when the interaction is stronger the energy 
levels of the machine become broadened, so even a sin- 
gle virtual qubit in the machine can couple to many vir- 
tual temperatures in the composite bath. Nevertheless, 
the principle remains the same; the external system is 
placed in thermal contact with the various virtual qubits 
in the bath, mediated by the virtual qubits of the ma- 
chine. 



V. REFRIGERATORS, HEAT PUMPS AND MORE 

As we already noted above, the particular function 
that the thermal machine provides depends upon the 
range of virtual temperatures of the composite bath 
which are coupled to the external system. If these are 
predominately lower than Ti and T2 (the actual temper- 
atures of the two baths) then the machine is functioning 
as a refrigerator. On the other hand, if the range of vir- 
tual temperatures are predominately higher than Ti and 
T2 then the machine is functioning as a heat pump. 

A more interesting situation is what happens when 
the range of virtual temperatures is predominantly neg- 
ative. This corresponds to virtual qubits having popula- 
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FIG. 3: Schematic diagram of smallest machine interacting 
with an isolated external qubit. An isolated external system 
- here a single qubit - is placed into thermal contact with the 
machine virtual qubit, which has matching energy-level spac- 
ing. This interaction is depicted by the wavy line. The net 
effect is that the external system is placed in thermal contact 
with virtual qubits in the composite bath at temperature T^, 
mediated by the machine virtual qubit. 




machine external 
virtual system 



/ 



FIG. 4: Schematic diagram of smallest machine interacting 
with an open external qubit. If the external system also inter- 
acts with its own environment at temperature T3 then it will 
behave as any system in contact with two thermal baths at two 
differing temperatures. 



tion inversions. The thermodynamic significance of this 
and related effects is discussed in the next section. 



VI. WORK AND HEAT ENGINES 

Work is one of the central concepts of physics. In 
comparison, population inversion is an interesting phe- 
nomenon, but until now usually associated mostly with 
laser physics. In this section however, we argue that in- 
version is a far more important thing - in fact we argue 
that 

Producing work and generating population in- 
version are one and the same thing. 

Work is generally associated with ordered movement. 
Therefore the first requirement in order to be able to talk 
about work is to have a system with many states, oth- 
erwise the whole idea of the system having an ordered 
movement i.e. going from one state to another in a sys- 
tematic way, makes no sense. Talking therefore of work 



when dealing with systems with a small, limited, num- 
ber of states, such as when considering the two qubits 
comprising the thermal machine presented in the previ- 
ous section, may seem impossible. To tackle this prob- 
lem many alternative ways of defining work have been 
given in the literature Il32ti35l , each with their own mer- 
its. Here we take the most direct route: just as Carnot 
said ''motive power (work) is the useful effect that a motor 
is capable of producing. This effect can always he likened to 
the elevation of a weight to a certain height." [36]. Hence 
we consider lifting a weight to be producing work, and 
any machine that lifts a weight to be a machine that pro- 
duces work. Following this idea we will argue that work 
is nothing other than producing population inversion. 

At a very elementary level, this definition is clearly 
consistent with the usual one. Indeed, consider a weight 
taken initially to be on the floor - this is the state of low- 
est energy and is therefore its ground state. Lifting the 
weight - that is, doing work on it - is nothing but com- 
pletely depopulating the ground state and populating 
an excited state, a state where the weight is higher up. 
Thus the final state of the weight is an inverted state. 
There is however far more here than meets the eye. 

To begin with we introduce the model of the smallest 
heat engine ||27l|3ll. It is identical to that of the refrig- 
erator, (i.e. we take 2 qubits, one connected to a bath 
at Ti and one connected to a bath at T2). We imagine 
that the engine delivers work by pulling up a weight. 
The weight is isolated from both baths. To simplify the 
situation we consider that the weight is pulled up very 
slowly, so that we can neglect the change of its kinetic 
energy and consider only the potential energy. Further- 
more, we suppose that the weight can be situated only 
at some discrete equidistant heights, so that the energy 
difference between them is the same. Hence, the weight 
has discrete energy eigenstates |n)^ with corresponding 
energy eigenvalues En = nE^, with £^ > 0. Alterna- 
tively, we can imagine that the engine delivers work by 
pumping energy into a harmonic oscillator; both situa- 
tions are formally almost equivalent (the harmonic os- 
cillator energies being limited from below). 

The free Hamiltonian of the system is thus 



Ho = Ei|l)/l| + E2|l),(l|+ «£^l«)>l (12) 



The energies are taken such that 

£2 - £1 = E^. (13) 

Given this constraint, the energy levels |0)Jl)2|n)^ and 
|1)J0)2 + are degenerate. The engine acts by mak- 
ing transitions between these degenerate states. 
The qubits and weight interact via the Hamiltonian 

00 

H/nf-g E |0)i|l)J»Li(l|2(0L(n+l| 

n=—oo 

+ |l)jO),|n+lL/0|,(lL(n| (14) 
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FIG. 5: Schematic diagram of the smallest heat engine, (a) The weight - an equi-spaced system, unbounded from above and 
below, interacts with the machine qubits via the interaction depicted by the orange and green arrows. This interaction induced 
transitions between the degenerate energy eigenstates |0)Jl)2|n)^^ ^ |1)J0)2 |n + l)^. The energies Ei and £2 ^re chosen such 
that the 'forward' transition in which the weight is Hfted is biased in favour of the 'backward' transition in which the weight is 
lowered. As such the weight is lifted on average, and hence this system produces work, (b) From the viewpoint of the virtual 
qubit, the biasing condition says that the virtual qubit must have a population inversion, or in other words must be at a negative 
temperature. 



The intuitive idea behind the design of this engine is 
to bias the transition 

|0)Jl)jM),„^|l)JO)Jn+lL (15) 

in which the weight is lifted in favour of the reverse 
transition in which the weight is lowered. This is ob- 
tained by coupling the two qubits to heat baths at differ- 
ent temperatures, T2 > Ti, chosen such that the proba- 
bility for the qubits to be initially in the state |0)Jl)2 is 
larger than the probability to be in the state |1)J0)2. 

In terms of the virtual qubit of the machine ^ what 
we want is to bias the transition 

|l)JnL^ |0)Jn + lL (16) 

in favour of the reverse transition. The condition for this 
to happen is simply that the probability to be in the state 
should be larger than the probability to be in the 
state |0)^. 

In other words, the condition for the machine to work as 
a heat engine is precisely that the machine virtual qubit has 
a population inversion and is therefore at a negative virtual 
temperature. 

A. Engine Model details 

So far our conclusions have been of a general nature; 
while they are enough to understand qualitatively the 
working of our heat engine it is illuminating to under- 
stand it also more quantitatively. To do this we need 
to say more about the way in which the two machine 
qubits interact with their respective heat baths. In Ap- 
pendix |B] we consider an explicit model for the thermali- 
sation oF the machine qubits, identical to that previously 
presented for the refrigerator [27J. This allows us to an- 
alytically solve for the time evolution of the system. We 
shall present here the important aspects of the solution. 



To begin with we must first introduce two new quan- 
tities relating to the virtual qubit. The first is the equilib- 
rium "bias'' of the virtual qubit, 

(z:^)=,(0|ti^T2|0),-,(1|ti^T2|1), 

ocl - e-^-/^- (17) 

which gives the difference in population of the ground 
and excited state of the virtual qubit at thermal equilib- 
rium (i.e. in the absence of interaction with the external 
system). Thus whenever there is more population in the 
virtual ground state then the bias is positive and so too 
is the virtual temperature. However, when we have a 
population inversion the bias and virtual temperature 
become negative. 

Second, we will need the equilibrium normalisation 
of the virtual qubit, 

«^)=,(0|Ti^T2|0),+,(1|Ti^T2|1), 

^ /l|2(0|Ti ^T2|1)J0)2 + ,(0|2(l|Ti ^T2|0)J1)2 

ocl + e"^^/^^ (18) 

which is simply the combined population in the virtual 
ground and excited states. Having introduced these two 
quantities we can now state the important features of the 
solution. To see that the weight is being raised in time, 
we must look at the average energ y of th e weight, 
Asymptotically we find that, from ( |B11| ) 

(E^) = -aE^{Zl^)t (19) 

2 

where oc = ^^f ^^2 is a positive constant which depends 
only upon the strength of the interaction Hamiltonian 
(g) and the strength of the thermalisation (p). We thus 
see that the energy of the weight is exactly propor- 
tional to the bias of the virtual qubit (Z^^), and thus in- 
creases in time whenever the virtual qubit is at a nega- 
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tive temperature, confirming our claim that the ''weight 
is lifted". 

Notably, we find that the weight is not only lifted in 
time, but that it spreads also. The expression obtained 
for the spread in the asymptotic limit is, from \B16) 

AEl = El{a{Nt^) - li{ZlY)t (20) 

where j6 = ^^p^i^^^2^p ^ is a second, model-dependent 
positive constant. 

Given the above, we see that the behaviour of the 
weight can be likened to that of a biased random walk: In 
this picture the virtual qubit plays the role of the coin, 
with the 'normalised bias' (Z^^)/ (N^^) playing the role 
of the (average) bias of the coin, and (N^^) modulating, 
along with g and p, the probability per unit time to make 
a transition. 



VII. STRENGTH OF WORK 

Classically when one thinks of work, one thinks of 
work full stop - there is no notion of there being differ- 
ent types of work. However, as we have seen above, all 
that is required from a heat engine in order to produce 
work is that its virtual qubit should be at a negative tem- 
perature. But different heat engines may have their vir- 
tual qubits at different negative temperature. For con- 
ciseness we will refer to this as to the "temperature at 
which the work is delivered" or the "temperature of the 
work". It is quite clear that work delivered at differ- 
ent temperatures must be in many physical ways differ- 
ent. In fact, intuitively, work delivered at negative tem- 
peratures with smaller absolute values - corresponding 
to larger inversions in the machine - is in some sense 
"stronger". 

It is important to emphasise that this "strength" of 
work is not equivalent to power, i.e. with the rate at 
which work is delivered; it is a completely different no- 
tion and has nothing to do with time scales. True, as 
we have seen in ([19| everything else being the same, the 
power of a machine depends on the temperature of the 
work - the smaller in absolute value the negative tem- 
perature, the larger the power. But power depends upon 
many other factors besides this temperature; critically, it 
depends on the coupling constants which are the ones 
that encode the time scales. On the other hand the vir- 
tual temperature at which work is delivered is a notion 
independent of any time considerations. 

The question that arises is how does the difference in 
virtual temperature of the work manifest itself physi- 
cally, apart from its influence on power? We do not yet 
know all of the ways in which it does, but there are two 
important scenarios which we shall demonstrate below. 



A. Producing inversions in systems bounded from above 

In the previous section we saw that the effect of a heat 
engine is to steadily increase the energy of the system on 
which it acts. This, of course, can only happen if the sys- 
tem's energy spectrum is unbounded from above. But 
suppose that the energy spectrum has an upper bound; 
what happens then? Intuitively the answer is immedi- 
ate: As we argued, all that a thermal machine does is to 
put the external system in simple thermal contact with 
the machine virtual qubit (which in turn is in thermal 
contact with a selected virtual qubit in the composite 
bath). If the external system is otherwise isolated it will 
tend to reach thermal equilibrium with the machine vir- 
tual qubit. If the machine virtual qubit is at negative 
temperature, this will be the final temperature of the 
external system. In a system with energy unbounded 
from above, such an equilibrium cannot be achieved - 
in the evolution towards equilibrium, the system will 
forever increase its energy. On the other hand, when 
the energy spectrum is bounded from above, the sys- 
tem will reach equilibrium, thermalising to an inverse 
Boltzmannian distribution at the negative temperature 
(see Appendix [A|. 

In other words, if the weight is lifted by a string pass- 
ing through a hole in a ceiling, the weight cannot raise 
more and will eventually reach a Boltzmannian distribu- 
tion, with the maximum probability at the ceiling, but 
having non-zero probability to be found at all smaller 
heights. 

Here we see one other aspect of the "strength" of 
work: The smaller in absolute value the negative tem- 
perature at which work is delivered, i.e. the "stronger" 
the work is, the closer to the ceiling the weight is 
pushed. 

The mechanism by which the inverted Boltzmannian 
distribution gets established is interesting by itself. To 
better see what happens consider an external system 
with many energy eigenstates. Suppose the system 
starts at an energy level far from the top. At the initial 
moment the top plays no role, and the engine simply 
pushes it up the energy ladder. Had this evolution been 
uniform, without spreading, the system would have 
simply ended "stuck at the ceiling", i.e. at the top en- 
ergy level. But as it climbs, the position of the weight 
also spreads - as stated previously, the evolution resem- 
bles closely that of a biased random walk. If initially the 
energy was perfectly well defined, during the evolution 
it will spread approximately as a gaussian packet. This 
is still different from the inverted Boltzmanian that is the 
final equilibrium state. However, as the system reaches 
the top energy level, it cannot climb in energy any fur- 
ther. At this moment, all it can do is to spread. As such 
it starts to extend backwards, toward lower energies, 
reaching eventually the inverted Boltzmanian equilib- 
rium state. 

Note that the spreading, that seemed to be more like a 
minor side effect when studying the evolution of a sys- 
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tern with infinite energy spectrum (or when far from the 
top in the case of a bounded energy one), plays in fact 
an essential role - the inverted Boltzmanian equilibrium 
could not be achieved without it. 

It is important to note also that the spreading is not re- 
stricted to quantum heat engines; classical thermal ma- 
chines lead to a spread in energy as well. Indeed, al- 
though usually one doesn't think of fluctuations in a 
thermal machine, they always exist. The thermal baths 
that drive the machines always have fluctuations and 
these lead to fluctuations in the evolution of the ma- 
chine. (A nice and easy example is given by the famous 
Feynman ratchet-and-pawl heat engine [37]. When the 
pawl is up the ratchet has some probability to move 
backwards.) 



B. Energy gain versus energy spreading 

Consider again a heat engine acting on system with 
an infinite number of energy levels and suppose that the 
system starts wit a well defined given initial energy. As 
we discussed before, it's evolution (in energy) is simi- 
lar to a biased random walk: as its average energy in- 
creases, it also spreads in energy. An interesting prob- 
lem is to compare the increase in average energy with 
the spreading. 

The increase in average energy is proportional to time, 
while the increase in spread is proportional to the square 
root of time (see equations ([19| and (|20|)). Thus for short 
times it is the spreading which dominates, so the sys- 
tem is likely to also lose energy, not only to gain it. The 
probability to find the system with lower energy there- 
fore increases, even though the average energy becomes 
larger. After a longer time however, the gain in the av- 
erage energy becomes larger than the spread. The in- 
teresting question is not after what time the average en- 
ergy becomes larger than the spread, but how much en- 
ergy must have been put into the system by this time. In 
other words, to measure the spreading versus the scale 
of average energy (as opposed to in time), which is an 
intrinsic property of the evolution. 

There is a "break even'' time f^^, such that 
{Ezv){ti)e) = ^^w{he)- ^6 denote the average en- 
ergy at this time (the "break even energy") by 
From ( p^ and j2Q| ) we find 

(21) 

The first term in the brackets depends only upon 
the virtual temperature, and not on the other de- 
tails of the machine. Indeed, it is easy to show that 
(]V^^)/(Z^^) = l/tanh(£^/2T^). The second term 
however is a model-dependent term, and thus the break 
even energy is a non-universal property of a thermal 
machine. 



Crucially however, when the virtual temperature ap- 
proaches minus infinity, i.e. the inversion of the ma- 
chine virtual qubit becomes vanishingly small, the sec- 
ond term becomes negligible and the break even energy 
becomes universal and infinite. This is another mani- 
festation of the strength of the work. Namely, when the 
strength is very weak, the break even energy becomes 
very large, and in the limit unobtainable. The machine 
therefore effectively stops working. As we will see in 
the next section, this is exactly the moment when the 
machine approaches the Carnot limit. 

That heat engines become "weak" in the Carnot limit, 
in that the power they deliver becomes vanishingly 
small, is a well known property. Here however, we 
presented a novel aspect, independent of time scales. 
We see that the Carnot limit is weak in that the weight 
must gain an infinite amount of energy before winning 
against its spread. Thus this is intrinsically a pathologi- 
cal point. 



VIII. UNIVERSALITY OF EFFICIENCY OF SMALLEST 
THERMAL MACHINES 

In the introduction we argued that the smallest ther- 
mal machines are the most fundamental thermal ma- 
chines, and because of this their study shows the inner 
workings of thermodynamics in the clearest way. Here 
we will show the insight they provide into the efficiency 
of thermal machines. 

The Carnot efficiency of a thermal machine is a uni- 
versal property of the machine - it depends only upon 
the temperature of the heat baths between which it 
works, but not upon the details of the construction of 
the machine or its interactions with the thermal bath. 
On the other hand, the efficiency of a thermal machine 
when operating away from the Carnot efficiency is a 
non-universal property, which does depend upon all of 
these details. In this section we will show however that 

the efficiency of the smallest machines is always 
universal - both at the Carnot limit and away 
from it. 

We shall focus here on the smallest refrigerator. Identi- 
cal conclusions follow for the heat pump and engine. 
The idealised thermal machine which the smallest re- 
frigerator is a specific instance of is depicted in FIG. [6] 
(a). The efficiency rj^^ for such a refrigerator is the ratio 
of the heat Q3 extracted from the 'cold' external bath at 
temperature T3 to the heat Qi provided to the machine 
by the 'hot' bath at temperature Ti, 

For the smallest machines it was shown previously 1381 
that the ratio of the heat exchanged with each bath must 
match perfectly the ratio of the energy-level spacing of 



9 



the corresponding qubits. 



Qi ■■ Qi ■■ Qs = El : E2 : E3 



(23) 



which arises due to the fact that there is only a single 
way in which the baths can exchange heat. Therefore 
the efficiency //^^ of the smallest refrigerator is given by 



q £3 



(24) 



Using the definition of T^, we can express this result not 
in terms of the energies Ei and £3, but instead in terms 
of the virtual and bath temperatures. We find that 



(25) 



where 



Vh — 1 



T2 



(26) 



is the standard Carnot efficiency for a refrigerator driven 
by baths at Ti and T2 and extracting heat from a bath at 
temperature T3 (see FIG. |6](a)). In section [IX| we will 
discuss this result with respect to approaching Carnot 
efficiency. 

In a similar fashion, we can derive the efficiencies t]^^ 

of the smallest heat pump and ?]^^ of the heat engine, 
given by 



T3 



where 



(27) 
(28) 



7^3 



- 1 - ^ 



(29) 



are the Carnot efficiency for a heat pump and heat en- 
gine driven by baths at Ti and T2, and delivering heat to 
a bath at T3, or work, respectively (see FIG.j6](b) - (c)). 

Thus we see that for the smallest machines the effi- 
ciency is universal - it depends only upon the tempera- 
tures of two external baths and the single virtual tem- 
perature in the composite bath which the machine se- 
lects and couples to. It does not depend on any other 
details, such as the coupling constants between the ma- 
chine and the bath, and the machine and the external 
system. 

The lack of universality of the efficiency of other ther- 
mal machines can now be better understood from this 
point of view. A general thermal machine will not select 



out a single virtual temperature, but rather a range of 
virtual temperatures, and couple an external system to 
all of them. The efficiency is therefore an average, which 
depends upon the relative coupling to each of these vir- 
tual temperatures, which is a machine-dependent and 
non-universal property. That is, it is precisely the abil- 
ity to 'perfectly filter' a single virtual temperature which 
gives the smallest machine its universal character. 



IX. REACHING CARNOT EFFICIENCY 

Previously in our study of the smallest thermal ma- 
chines we studied the question of whether they can ap- 
proach the Carnot efficiency. Naively one may think that 
imposing restrictions upon the size of a machine may 
necessarily constrain it so that it cannot, even in prin- 
ciple, operate anywhere close to the Carnot efficiency. 
This however is not the case; as shown previously 1 28 J 
the smallest machines can approach the Carnot effi- 
ciency. What we will show here is a surprising twist 

Essentially only the smallest machines can ap- 
proach Carnot efficiency. 

(a more detailed and qualified statement shall be pre- 
sented below). 

In the previous section we showed that the efficiency 
of the smallest machines can be written as the Carnot 
efficiency plus a correction term, equations ( |25| ) - ( |28| ). 
To approach the Carnot limit this correction term must 
vanish. Hence the machine must be fine-tuned to cou- 
ple to the specific virtual temperature that makes this 
term vanish. For a machine constructed by two physical 
qubits, this may be achieved easily by tuning the energy 
level spacing Ei and £2 appropriately. 

Let us consider now more general thermal machines. 
As stated previously, unlike the smallest thermal ma- 
chines, which filter out a single virtual temperature from 
the composite bath, generally a thermal machine will 
filter out a range of virtual temperatures. The above 
analysis shows that if a machine filters out a tempera- 
ture away from the precise one required, then it will not 
operate close to the Carnot limit. Thus in order to ap- 
proach Carnot efficiency general thermal machines must he 
engineered to decouple from the un-wanted virtual tempera- 
tures, i.e. the machine must filter out a range of virtual 
temperatures which approach the single desired tem- 
perature in the limit. 

The only way to achieve this de-coupling is for the 
machine to be engineered in such a way that it effec- 
tively implements a smallest thermal machine. Indeed, 
consider a single pair of energy levels at spacing £3 in 
the external system which the machine is acting upon. 
Given fixed external bath temperatures T\ and T2 we 
noted previously that there is only a single choice of 
spacings E\ and £2 that both couple to £3, (i.e. E^ — Ex — 
£3) and also produce a machine virtual qubit at temper- 
ature Ty. On the other hand, if the parts of the machine 
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in contact with each bath have other spacings E[ and E2 
which can couple to £3 (i.e. — E[ = £3) then they 
lead to other virtual temperatures different from the 
required T^. Coupling the machine to these other virtual 
temperatures will spoil its efficiency. Thus we conclude 
that any thermal machine which filters a single virtual tem- 
perature must he coupled by only a single spacing to each bath 
and is hence exactly the smallest thermal machine or a collec- 
tion thereof. 

To summarise, quantum mechanics offers for free sys- 
tems with discrete energy levels. All we need to do is en- 
sure the spacing is the desired one. On the other hand, 
classical systems have essentially a continuous energy 
spectrum, so we cannot avoid them having undesired 
energy level spacings. It is nevertheless possible to ap- 
proach Carnot efficiency if we engineer the machine in 
such a way that the undesirable spacings that couple to 
the system are all close to the desired virtual tempera- 
ture. Hence classical machines that approach the Carnot 
limit do so by effectively removing energy levels from 
their spectrum, and becoming essentially identical to the 
smallest machines. ll4H 

Finally, we would like to comment on the meaning 
of the Carnot limit, and the way in which reversibility 
is achieved. As it is well known, whenever two bod- 
ies at two different temperatures are put in contact, irre- 
versibility occurs. The only way to achieve reversibility 
is to ensure that all systems that are in contact have the 
same temperature. This is precisely what happens in a 
machine which works in the reversible regime. As we 
can see from ( |25| ) and ( |27| , for the refrigerator and heat 
pump the correction term vanishes and the Carnot limit 
is achieved precisely when the bodies which are placed 
into thermal contact - the external system, the machine 
virtual qubit, and the selected virtual qubits of the com- 
posite bath - all approach the same temperature. Put 
differently, one may be tempted to think that we have 
a situation in which three bodies at different tempera- 
tures are in contact with each other and this may nec- 
essarily lead to irreversibility. The difficulties that may 
arise from such a situation were exposed nicely by Par- 
rondo 1391 . We see however that reversibility can occur 
even in such a situation precisely because the relevant 
systems which interact in this situation are both at the 
same temperature. 



X. ENTROPY 

Analysing the flow of entropy in thermal machines is 
one of the most basic aspects of thermodynamics. Here 
we will do so for the smallest thermal machines, where 
we find that it is highly illuminating. In particular we 
will see that entropically as well the virtual qubit be- 
haves as a real qubit at the virtual temperature, rein- 
forcing our argument that the virtual temperature acts 
as a real temperature in all respects. Let us focus on 
the smallest refrigerator, depicted in idealised form in 




QiNp Si 



Qi fsi 






FIG. 6: Schematic diagram of idealised thermal machines. 

Schematic diagram, showing the flow of heat and entropy for 
(a) refrigerator (b) heat pump (c) heat engine. 



FIG. [6] As noted previously, the smallest machines have 
the property that the ratio of heat currents match the ra- 
tio of energy levels, (|23j. Using this fact, in conjunction 
with the definition of T^, ([9|, we obtain 



Q2 _ Qi _ Q3 
T2 Ti 



(30) 



which is clearly an entropy equality involving the en- 
tropies Si = Qi/Ti, S2 = Qil^i, flowing out of bath 1 
and into bath 2, and the entropic quantity = Q3/T'y- 
The change in total entropy of the two baths and ex- 
ternal system for the refrigerator is given by 



(31) 



where S3 = Q3 / is the entropy flowing out of bath 3. 
Using ( (30| this can be re-written as 



ASf,= 



Q3 _ Q3 



(32) 



By similar analysis, for the heat pump and heat engine 
we find 



AC _ Q3 Q3 

^^hp — ^ 



AShe = 



Q3 



(33) 



The meaning of these results (for the refrigerator and 
heat pump) are as follows: We have two bodies in ther- 
mal contact - the external system at temperature T3, and 
the virtual qubits in the composite bath at temperature 
Ty. The bodies exchange the quantity of heat Q3. Now 
we see that the meaning of the entropic quantity Q3/Ty 
is the entropy of the virtual qubits, and the relation ( [32] ) 
is the standard entropy flow for two bodies in thermal 
contact. 

Similar analysis of a machine putting heat into a body 
at < shows how the strength of work, as charac- 
terised by the smallness of the absolute value of T^, is 
related to the amount of entropy that is removed from 
the body while delivering a fixed amount of heat Q3: 
the smaller the absolute value of the larger the en- 
tropy removed from the body, with the amount tending 
to infinity as ^ from below. 
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Incidentally, this shows again that the Carnot limit is 
approached exactly when the two systems in thermal 
contact approach the same temperature. 

XI. GENUINE THERMAL MACHINES 

We would like to define now an important concept, 
that oi genuine refrigerators and heat pumps. Your refriger- 
ator at home uses a supply of work to transfer heat from 
a cold source to a hot source. Furthermore, from a theo- 
retical point of view, the way refrigerators are generally 
presented is that of "an engine running in reverse'', us- 
ing a supply or work to move heat against the gradient. 
The point we want to make now is that work is not neces- 
sary. 

On the one hand it is well known that an external 
source of work is not required to produce a refrigerator. 
However, by itself this fact is not surprising - we can 
always replace the source of work with an engine that 
produces it, using two external baths. Indeed, we know 
that the work used to power our refrigerator at home 
ultimately comes from a power station. There are also 
examples of refrigerators in which there is no apparent 
place in which work is used - adsorption refrigerators 
for example have no engines and pumps - but this by 
itself doesn't mean that work is not present in some hid- 
den form (such as while gases are expanding, etc.) The 
concepts of virtual qubits and virtual temperatures offer 
for the first time the tools to address this question. 

The crucial point we want to make is that it is not nec- 
essary to generate work - even internally in the machine 
- to create a refrigerator. The refrigerators we presented 
in the previous section function by putting an external 
system in contact with virtual qubits whose temperature 
is lower than Ti and T2, the actually temperatures of the 
two baths. At no point is there ever a population in- 
version involved - in other words at no point is work 
produced or used. 

In this sense, this is a genuine refrigerator. Ordinary 
refrigerators (and heat pumps) that use work are now 
seen to be "wasteful". This is most evident if we talk 
about the resources that a composite thermal bath pro- 
vide us with. It is clear that work is a resource provided 
by the composite bath, and we can indeed use this to 
produce a refrigerator. However, we see that the com- 
posite bath provides other resources as well, namely 
energy at any temperature, not merely work. Energy at 
a cold temperature is a resource which can be used to 
achieve cooling directly, while work needs to be some- 
how converted. Thus ordinary refrigerators make non- 



necessary use of a resource (work) that is more powerful 
than what is needed. Genuine thermal machines use the 
minimal possible resources, and are thus the "purest" 
thermal machines. 

Furthermore, it is essential to make the difference 
between genuine and non-genuine thermal machines 
when we want to know what exactly do they do to the 
thermal baths that drive them. If we take the baths not 
to be strictly infinite in size, then thermal machines nec- 
essarily degrade them, in particular reducing their free 
energy. Genuine refrigerators couple to different virtual 
qubits than refrigerators which use work as an interme- 
diate effect. Hence they clearly affect the bath in a differ- 
ent way. It would be interesting to study further exactly 
these differences. 

Finally, it is illuminating to see, for the case of the 
smallest machines, how the functionality is changed as 
we vary the bath parameters. In FIG. [7| we hold fixed 
the design of the machine (£1 and E2), and the bath tem- 
perature Ti, and plot how the virtual temperature (and 
therefore function) changes as we vary the remaining 
bath temperature T2. It is interesting to note that as T2 
approaches from above then 00. This engine 

is the reversible (Carnot) engine, and is seen to be at the 
transition between heat pump and heat engine; in this 
sense it is again seen to be the 'weak'. 



















refrigerator 


heat 1 / 


heat 






pump I / 


engine 





FIG. 7: Graph of virtual temperature against bath tempera- 
ture We hold fixed the local energy-level spacings £1 and £2 
as well as the bath temperature T\, and allow the bath temper- 
ature T2 to vary. When T2 < T^, the virtual temperature 
becomes smaller than either environmental temperature and 
the machine is a refrigerator. For < T2 < |^Ti, T^:; becomes 
larger than T\ and T2 and hence the machine is a heat pump. 
Finally, for |^ < T2, becomes negative and the machine 
functions as a heat engine. 
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Appendix A: Thermal machines acting on finite 
dimensional systems 

In this appendix v^e v^iW demonstrate explicitly our 
claim, that external systems put into thermal contact 
w^ith the machine virtual qubit thermalise to its virtual 
temperature, for the case of finite dimensional systems. 
Crucially this holds for both positive and negative tem- 
peratures. The temperature which the external system 
thermalises to classifies the behaviour of the machine: If 
it is colder than either environmental temperature then 
the machine is a refrigerator; if it is hotter then it is a 
heat pump. Finally, if the temperature is negative, then 
it is a heat engine. 

To see explicitly this thermalisation, consider initially 
that the external system is itself completely isolated ex- 
cept for the interaction with the thermal machine. That 
is, we consider the external system to only be in contact 
with the thermal machine, not with any other external 
system, so it is only the machine that determines its be- 
haviour. 

In the weak coupling limit, the dynamics of the sys- 
tem are accurately described by a master equation. The 
equation will generically take the form 

1^ = -f [Ho + Hi„t,p] + Pi (|0) + V2{p) (Al) 

where 

Ho - Ei|l)/l| + E2|l),(l| + £311)3(11 (A2) 

Hint=g(\o\\i\\ouiwm 

+ |1)J0), 11)3/01,(113(01) (A3) 

with E3 = E2 — El/ and Vi (p) and ^2(p) are the dissipa- 
tive dynamics acting on qubit 1 and 2 respectively. Such 
an equation provides a consistent description of the dy- 
namics of the system when both the strength of the in- 
teraction Hamiltonian and the dissipative dynamics are 
weak 1401 . 

Without specifying a specific model for the dissipative 
dynamics, all such models must satisfy 

r>i(Ti ^ (723) = ^2(^2 ^ 0-13) = (A4) 

where 

Ti = ^e-^^'^' Zi = tr(e-^^/^0 (A5) 

are the thermal equilibrium state and partition function 
of qubit z, and the states ^23 and ^13 are arbitrary density 
matrices. In other words, a necessary requirement for 
any weakly coupled and weakly interacting system is 
that the thermal equilibrium state of each qubit is the 
stationary state of its dissipative dynamics. 
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To find the stationary solution to ^Al\ it suffices to 



note that 



Ml ^ El ^ El Hn H- — Q f ^ — ^ — 

Ti ^ T2 ^ T3 ' ^0 ^ J^int^ — 6 Ti To,) 

(|0)J1),|0)3/1|,(0|3(1| - |1)J0)J1)3/0|,(1|3(0|) (A6) 



Thus, by recalling the definition of the virtual temper- 
ature given in the main text, equation it is evident 
th at it is exactly when T3 = that the right-hand-side 
of \A6) vanishes. In this instance the operator appearing 
on the left of the commutator commutes with the total 
(free + interaction) Hamiltonian of of the system. If we 
thus define the state 



(A7) 



which is simply the thermal equilibrium state of qubit 3 
at temperature T^, then it follows immediately that the 
thermal product state 



Ti (g) T2 (8) 



111 



Hi 



Hi 
T2 



H3 



also commutes with the total Hamiltonian - showing 
that it is stationary with respect to the unitary dynam- 
ics. This state is however also seen to satisfy conditions 
( |A4| , showing that it is stationary with respect to the dis- 
sipative dynamics. Combining these two facts we thus 
conclude that it is the stationary solution of the master 
equation ( |A1| ). 

Thus the external qubit will approach a thermal state 
at the temperature of the machine virtual, qubit showing 
that it indeed thermalises to the temperature of the virtual 
qubit, as stated. 

It is important to realise that this result can be ex- 
tended beyond the thermalisation of qubits. Indeed, 
qubits have the special property that every diagonal 
density matrix (in the energy eigenbasis) can have a tem- 
perature associated to it. Thus to show that the exter- 
nal system is indeed thermalising to the temperature of 
the virtual qubit we must consider more general sys- 
tems. Let us consider therefore an N-level system which 
has equidistance spacing between its energy eigenstates. 
That is, we consider a system with energy eigenstates 
|0), . . . , |N — 1), and free Hamiltonian 



N-l 



(A9) 



n=0 



where we maintain the condition that E3 = E2 — Ei. 
This system is taken to interact with the two machine 
qubits via the interaction Hamiltonian 



«!2-^Elo)Ji).K(i|.(oi3(« + i| 

+ |l)j0),|n + l)3/0|,(l|3(n| (AlO) 



N-2 




FIG. 8: Schematic diagram of smallest machine interacting 
with an isolated external system. An isolated external system 
- here an N level equi-spaced system - is placed into thermal 
contact (wavey lines) with the machine virtual qubit, which 
has matching energy-level spacing. The net effect is that the 
external system is placed in thermal contact with virtual qubits 
in the composite bath at temperature T^, mediated by the ma- 
chine virtual qubit, and reaches a Boltzmannian at tempera- 
ture Ty. This holds independent of whether is positive, or 
negative, in which case the state is an 'inverted Boltzmannian'. 



which can be seen as a sum of terms, all of the form ( |A3| ). 
Using an identical proof method as used for qubits, or 
by direct substitution, it is straightforward to show that 
the thermal product state 



Ti (g) T2 I 



rim 



1 



exp 



Hi_^ 



(All) 

is the stationary solution to the corresponding master 
equation, where 



m _ 



(A12) 



In other words, all equi-spaced systems, when placed into 
therm al contact with the virtual qubit (via the interaction 
\A1Q)) ) thermalise to the virtual temperature. Crucially, 
when the virtual temperature is negative, the final state 
is an inverse Boltzmannian, with population decreasing 
exponentially with decreasing energy as we move away 
from the most excited state |N — 1). 

If instead of considering the external system to be iso- 
lated we place it in contact with a third thermal bath at 
some fixed temperature T3 then we would expect to find 
a non-universal steady-state solution to \A1) which will 
depend upon the explicit form the of dissipative terms 

In [27J, for the case where the external object is a qubit 
(to be cooled), a specific model was employed and it was 
shown that with three thermal baths the external object 
reaches a stationary state which is cooler than its envi- 
ronment. From the current perspective the situation is 
simple to understand; the object is in thermal contact 
with both the virtual qubit and also its own environ- 
ment. We thus expect it to reach a steady state with 
temperature between these two temperatures, with the 
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precise temperature depending on the relative 'magni- 
tude' of the two thermalising effects. Therefore as long 
as Ty < intuitively we expect the object to be cooled, 
which is exactly what was found for the model pre- 
sented in lEZl. 



Appendix B: Thermal machines acting on infinite 
dimensional systems 

In this appendix we will present the detailed calcula- 
tions of a specific model which shows that when we take 
the external system to be unbounded in energy, and en- 
gineer the smallest thermal machine to have the virtual 
qubit at a negative virtual temperature, then the system 
is continually excited, and thus captures the classical no- 
tion of a heat engine producing work, as claimed in the 
main text. 

As explained earlier, to do so we place the machine 
virtual qubit in thermal contact with a "weight" - a 
system consisting of equally-spaced energy levels, un- 
bounded from above and below. Namely we have the 
free and interaction Hamiltonians 



Hr 



= £i|l)/l|+£2|l),(l| 



nE^\nl{n\ (Bl) 



Hrnt=g E I0)i|l)2l^li(l|2(0|> + l| 

+ |l)JO)Jn + l),„/0|,(l|>| (B2) 

where Ei and E2 are chosen such that £2 ~ = To 
complete the model we must give an explicit form for 
the dissipative dynamics of the machine, i.e. a model for 
the thermalisation of each machine qubit with its bath. 
We will use the same model used in our previous work 
on the smallest refrigerator [27], namely 



(B3) 



This simple model of thermalisation clearly satisfies the 
conditions ( |A4| and and generates an exponential decay 
of each qubit back to its thermal state T/, with decay con- 
stant p. The explicit master equation which governs the 
dynamics of the system therefore takes the form 



dp 



-i[Ho + Hi„t,p] 



! = 1 



(B4) 



We expect that this equation accurately models the be- 
haviour of the system as long as g and p are chosen such 
that grP <^ so that we are in the weak-dissipation 
weak-interaction regime. Here it is consistent both to 
use the free Hamiltonian to define the thermal states, 
and to neglect additional contributions to the dissipa- 
tive dynamics originating from the interaction between 



the qubits. Note that here we make the simplifying as- 
sumption that the two baths interact with the same cou- 
pling strength p. The analysis presented below can also 
be carried out in the case of unequal coupling strengths 
pi ^ p2, however it substantially complicates the alge- 
bra without changing the results obtained. 



1. Raising of the weight 



By first multiplying ( [B4| ) by the free Hamiltonian of 
the weight H^, and then taking the trace, we find an 
expression for the rate of change of the average energy 
of the weight, given by 

j^{E^) = ^tr(H^p) = -igE^A (B5) 

where 

A = Ei(0l2(l|>IP|l)il0)2l" + lL 

n 

-/l|,(0|> + l|p|0)Jl)JnL (B6) 

By introducing the three new operators on qubits 1 and 
2, 

Z, = |1)J0),/1|,(0| - |0)J1),/0|,(1| 
Z, = |0)J0),/0|,(0|-|1)J1),/1|,(1| (B7) 
= |1)J0),/1|,(0| + |0)J1),/0|,(1| 

and denoting (Z^) = tr{^pZv), (Zy) = tri^pZ^) and 
(Nj,) = tri^pNv) then the following set of equations can 
be obtained in a straightforward manner from | |B4| | 



dt 



-2pA - lig{Zy) 



d_ 
dt ' 
d_ 
~dt 
d_ 
di 



(Z,) = -2fgA-p((Z,)-(z:^)) (B8) 

:(N.) - p(l-2(N,) + (Z^^)(Z,)-(Z;''')(Z;)) 

(z,) = -p((z;)-(z;^^)) 

where 



tr(ZyTi I 



(B9) 



are the thermal equilibrium expectation values of the 
operators Z^ and Z^. We see that the evolution of (Z^) 
is completely independent of the evolution of the other 
variables - a situation which only occurs when the two 
baths have a common coupling strength p. Furthermore, 
although the evolution of (Nj;) depends on all the other 
variables, we see that no other variable depends directly 
upon it. This system of equations can easily be solved 
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for the asymptotic (stationary) behaviour. Denoting the 
variables in this limit with a superscript 'S' we find that 



















(BIO) 









There are a number of notal^le features of this solution. 
The first is to note that (Z^ ) is the 'bias' of the "anti- 
virtual'' qubit - the virtual qubit formed from the states 
1 00) and |11), and that this bias is left unaltered by the 
coupling of the virtual qubit to the weight. Second, we 
see that the normalisation of the virtual qubit is shifted 
from its value at thermal equilibrium and that its sta- 
tionary value is only ever smaller than (or equal to) this 
value. Finally, we see that the bias of the virtual qubit is 
scaled by the coupling, and always decreases in magni- 
tude. 

Finally, but substitution of the solution \B10) into ( [B5| ), 
we obtain for the asymptotic rate of change of the aver- 
age energy 



(Bll) 



2. Spreading of the weight 

A second quantity of interest is which gives in- 
formation about the spreading of the average energy of 
the weight in time. From ( |B4( j it follows that the equa- 
tion governing the evolution of this quantity is 



dt 



(B12) 



where we have introduced the quantity A^, defined by 

00 

An= E (2n + l)UO|,(l|>|p|l)JO),|n + l)^ 

n=—oo 

-/l|,(OL(n + l|p|0)Jl),|n)^) (B13) 

Introducing a second new quantity E^, related in an 
analogous manner to (Z^) as An is to A, 

00 

X] (2n + l)UlL(0L(n + l|p|l)j0)Jn + l)„ 

n=— 00 

- mMJ.n\pmi\\nl) (B14) 

Then it is possible to obtain the following set of coupled 
differential equations 

^-An = -IpAn-ligLn (B15) 



dt 
d_ 
dt 



which are valid in the asymptotic limit, since to derive 
the second equation the asymptotic solution ( |B10| ) has 
been used. These equations can easily be solved using 
standard techniques and the asymptotic solution for the 
expected-squared-energy is given by 



-(£2' 



g2p 



2g2(g2+2p2) 
(2^2 + p2)2 



3. Heat transfers 



^^^)') (B16) 



Finally, if we define further the two quantities Fi and 
F2, which are the instantaneous ground state probabili- 
ties for qubits 1 and 2 respectively, 

ri = Ei(oi2(ol>lp|o>ilo>2l"L 

n 

+ /0|,(l|>|p|0)Jl)Jn),„ (B17) 
T2 = EPU0l{n\p\0\\0\\nl 

n 

+ ,{imMp\^\mni (BI8) 

then these quantities obey the coupled set of equations 

d 



dt 
d_ 
dt 



Fi = +zgA + p(F'i^-Fi 
F2 = -zgA + p(F^^-F2 



(B19) 



where F^^ = (1 + is the equilibrium ground 

state population of each qubit. Making use of (b10| ) it 
can therefore be seen that asymptotically these popula- 
tions reach the values 



yS _ Y^q g^ (rycq^ 

^2-^2 ~ 2g2+p2 [^v ) 



(B20) 



The rate at which heat flows between the qubits and 
their environments is given by the change in energy of 
each qubit due to the interaction with the baths. Given 
the master equation ( [B4] ), the asymptotic heat currents 
are therefore 



|q,- = tr(H,P,(p)) = ptr^Hiin-pf)) 

2 



(B21) 



